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Abstract: On a deglaciated terrain, glacial gravel is the primary component of the natural habitat for 
vascular plant colonization and succession. Knowledge regarding the role of glacial gravel in vascular plant 
growth, however, remains limited. In this study, an unmanned aerial vehicle (UAV) was used to investigate 
plant family composition, species richness, fractional vegetation cover (FVC), and gravel cover (GC) along 
elevational gradients on the three glacier forelands (Kekesayi, Jiangmanjiaer, and Koxkar Baxi) of the 
Third Pole (including the eastern Pamir Plateau and western Tianshan Mountains) in China. We then 
analyzed the spatial characteristics of vascular plants followed by exploring the effect of glacial gravel on 
vascular plants. Findings indicated that FVC on these glacier forelands generally decreased as the elevation 
increased or distance from the current glacier terminus decreased. The shady slope (Kekesayi) was more 
vegetated in comparison to the sunny slope (Jiangmanjiaer) at the glacier basin scale, and the warm and 
humid deglaciated terrain (Koxkar Baxi) had the highest FVC at the regional scale. Plant family 
composition and species richness on the glacier forelands decreased with rising elevation, with the 
exception of those on the Jiangmanjiaer glacier foreland. The relationships between FVC and GC 
presented negative correlations; particularly, they exhibited variations in power functions on the Kekesayi 
and Jiangmanjiaer glacier forelands of the eastern Pamir Plateau and a linear function on the Koxkar Baxi 


glacier foreland of the western Tianshan Mountains. Glacial gravel was found to be conducive to 
vegetation colonization and development in the early succession stage up until vascular plants adapted to 
the cold and arid climatic condition, whereas it is unfavorable to the expansion of vascular plants in the 
later succession stage. These findings suggested that the spatial difference of plant characteristics had 
close connections with regional climatic and topographic conditions, as well as glacial gravel distribution. 
In addition, we concluded that aerial photographs can be an asset for studying the functions of 
micro-environment in vegetation colonization as well as succession on the glacier forelands. 
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1 Introduction 


There are a large number of mountain glaciers in China. According to the second Chinese Glacier 
Inventory (CGI) (Liu et al., 2015) and Randolph Glacier Inventory (RGI) 6.0 (RGI Consortium, 
2017), the country's glaciers cover 5.18x10* km? or 7.3% of the world's total glacial area. This is 
the highest percentage with the exception of the Arctic and Antarctic (70.58x10* km°). Due to 
global warming, glaciers worldwide have typically been retreating since the end of the Little Ice 
Age (approximately 1850) (Leclercq et al., 2011; IPCC, 2012), and China is no exception (Li et 
al., 1998; Liu et al., 2003; Shangguan et al., 2006; Zhang et al., 2016; Li et al., 2019). Uncovered 
surfaces after glacier retreat offer excellent opportunities to examine the processes of primary 
vegetation succession in the field. As a result, a number of research findings on succession 
dynamics along with factors impacting vegetation development have been published, which are 
primarily implemented in High Arctic, European mountainous areas, and North America (Raffl et 
al., 2006; He and Tang, 2008; Burga et al., 2010; Wietrzyk et al., 2018; Wei et al., 2021). In 
China, the prime focus of the studies on the spatial distribution and succession of vegetation is on 
several glacier forelands, which consist of Yulong Snow Mountain, Tianshan Mountains, and 
Gongga Mountain (Li and Xiong, 1995; Chang et al., 2014; Wei et al., 2021). These works have 
provided valuable vegetation information on the glacier forelands of the Third Pole region in 
China. However, in comparison to those in the aforementioned regions, there is a scarcity of 
related studies in this field (Third Pole region). 

All stages of plant community development on the glacier forelands are dominated by both 
biotic and abiotic factors, whereas in the early stage of vegetation succession, abiotic control 
plays a more crucial role (Houle, 1997; Jumpponen et al., 1999). Glacial gravel, which comes 
either supraglacially from nunataks and valley sides or from erosion of the subglacial bed 
(Boulton, 1978), is quite prevalent on the glacier forelands. Scientists have taken a keen interest 
in it since vegetation tends to grow in surrounding glacial gravel. According to the previous 
publications (Rooney, 1997; Mong and Vetaas, 2006), glacial gravel plays a controversial role in 
vegetation colonization as well as the development of the glacier forelands. In particular, some 
studies have pointed out that glacial gravel on the glacier forelands prevents disturbances from 
herbivores (Rooney, 1997) and provides relatively safe sites for the early colonization of several 
pioneer plants (Stécklin and Bäumler, 1996; Jumpponen et al., 1999; Niederfriniger Schlag and 
Erschbamer, 2000). While on the other hand, other studies have argued that the potentially 
protective glacial gravel at a microsite becomes an obstacle for seedling establishment because 
the available space for germination is occupied (Mong and Vetaas, 2006). Therefore, there is an 
urgent need to further explore the spatial changes of vegetation characteristics as well as the 
impact of glacial gravel on vegetation colonization and development. In order to comprehensively 
understand this scientific question, our investigation had to be expanded to several glacier 
forelands with various regional typographic and climatic conditions. 

Glacier forelands are characterized by complex terrains, unique ecosystems, and harsh 
environments, making the field investigation of vegetation dynamics more challenging. Up to 
now, traditional ground sampling has been a valuable approach and has been extensively applied 
in field surveys (Niederfriniger Schlag and Erschbamer, 2000; Mong and Vetaas, 2006; Raffl et 
al., 2006; Burga et al., 2010; Wietrzyk et al., 2018), but this methodology is time-consuming and 
labor-intensive, and requires trampling fragile habitats (Lee and Yeh, 2009; Martinez-Lopez et al., 
2014). Additionally, the free satellite image is difficult to utilize for the assessment of vegetation 
characteristics on the glacier forelands owing to its low resolution. In contrast, the unmanned 
aerial vehicle (UAV) provides an ideal solution to this problem due to its natural advantages, such 
as increased field efficiency, effective avoidance of human disturbance, the automatic acquisition 
of high-resolution images at the quadrat scale, and the ability to repeatedly observe the same area 
after setting the flight line (Sun et al., 2018). Furthermore, the UAV can easily acquire vegetation 
data in inaccessible areas. These advantages make the UAV technique suitable for vegetation 
monitoring in a variety of topographic conditions including glacier forelands. 
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Fractional vegetation cover (FVC), plant family composition, species richness, dominant 
species, and gravel cover (GC; the ratio of the vertical projected area of glacial gravel to the total 
ground area) along chronosequences were investigated by integrating ground-based sampling and 
aerial photography on the three glacier forelands of China's eastern Pamir Plateau and western 
Tianshan Mountains in this study. Following that, we analyzed the variation characteristics of 
these vegetation indices during primary succession along with the relationship between GC and 
FVC. This study aims to characterize the spatial variations of vegetation characteristics (FVC, 
plant family composition, and species richness), examine the spatial heterogeneity of vegetation 
under various regional typographic and climatic conditions, and assess the influence of glacial 
gravel on vegetation community development. The findings of this study will assist investigators 
of primary vegetation succession and biodiversity, and aid in regional ecological conservation and 
sustainability in the Third Pole region. 


2 Materials and methods 


2.1 Study area 


The study area is situated in the Third Pole region and spans across the eastern Pamir Plateau and 
western Tianshan Mountains in Xinjiang Uygur Autonomous Region, China. A total of three 
glacier forelands namely Kekesayi, Jiangmanjiaer, and Koxkar Baxi were surveyed. Within each 
glacier foreland, the sampling sites were established along elevational gradients (Fig. 1). The sites 
in the eastern Pamir Plateau were distributed in the retreat areas of modern and paleo glaciers 
(elevation ranges of 3843-4004 m a.s.l. for the Kekesayi glacier and 3631-4241 m a.s.l. for the 
Jiangmanjiaer glacier), whereas those in the western Tianshan Mountains were primarily chosen 
in the deglaciated terrain of the modern glacier (elevation range of 2979-3099 m a.s.l. for the 
Koxkar Baxi glacier). 

The Kekesayi glacier (code number of 5Y663D0087 from the Chinese Glacier Inventory) and 
Jiangmanjiaer glacier (code number of 5Y663D0004 from the Chinese Glacier Inventory), which 
are typical valley glaciers, are situated in the same glacier basin (Fig. 1b—d) under the ice caps of 
the Muztagata Peak and Kongur Peak, respectively. With an approximate length of 18.0 km and 
an extent of 86.50 km?, the debris-covered Kekesayi glacier extends from the east to the north and 
is by far the largest glacier of the Muztagata Peak (Shangguan et al., 2006; Seong et al., 2009a, b). 
The Jiangmanjiaer glacier, which is 14.8 km long and extends over 45.08 km?, is the longest 
glacier on the southwestern slope of the Kongur Peak (Wang et al., 2011). The field investigation 
of the moraine ridge and previous publications on the quaternary glacier landform revealed that 
the entire Muztagata Peak and Kongur Peak were covered by huge glaciers in the historical 
period, and glaciers even extended into the Kara Kol Lake during the glacial maximum of these 
areas (Cui, 1960). However, based on remote sensing data, their frontal positions have had almost 
no visual changes in the last four decades (Holzer et al., 2015). The Muztagata Peak and Kongur 
Peak are mainly controlled by the upper-level westerly circulation and local circulation, and 
moisture transference is difficult here because of the "blocking effect" of surrounding high 
mountains (Yu et al., 2006). As a result, the two peaks exhibit a cold and semi-arid continental 
climate. The Taxkorgan meteorological station, situated approximately 50.0 km south of the 
Muztagata Peak, is the sole station above 3000 m a.s.l. on the eastern Pamir Plateau. From 1957 
to 2010, the mean annual temperature at this station was measured as 3.4°C, the mean summer 
temperature as 15.1°C (June—August), and the mean annual precipitation as 70.2 mm (Yan et al., 
2013; Yang et al., 2014). The distribution of vegetation types along elevational gradients is as 
follows: meadow (below 3500 m a.s.l.), desert vegetation (3500-4000 m a.s.1.), alpine grassland 
vegetation (above 4000 m a.s.l.), and alpine sparse vegetation (above 4500 m a.s.1.); the floristic 
composition consists of Ceratoides, Graminoids, Artemisia, Stipa, etc. (Wang et al., 2016). The 
vegetation characteristics of the two glacier forelands, nevertheless, were obviously different (Fig. 


2a-f). 
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Legend 
Sampling site — Contour line C3 Glacier WE Glacier basin EE Other region 


Fig. 1 Overview of the study area (a) and distribution of sampling sites in the Kekesayi glacier (b) and 
Jiangmanjiaer glacier (c) of the eastern Pamir Plateau (d) and in the Koxkar Baxi glacier of the western Tianshan 
Mountains (e) 


The Koxkar Baxi glacier (code number of 5Y674A0005 from the Chinese Glacier Inventory), 
situated on the southern slope of the Tuomuer-Khan Tengri Mountain in the western Tianshan 
Mountains (Fig. 1), is a typical continental glacier with the length of about 26.0 km, an area of 
83.60 km’, and elevations of 3020-6342 m a.s.l. (Zhang et al., 2007). It is a glacier covered in 
debris, and since the 1980s, its terminus has been retreating at an average rate of 0.5-1.5 m/a (Xie 
et al., 2007). The mid-latitude westerlies that originate from the Atlantic Ocean have the most 
impact on the region. The Koxkar Baxi glacier presents a sub-humid climate with the mean 
annual precipitation of 669.4 mm at 3000 m a.s.l. (Li et al., 2012). The mean summer temperature 
(June—August) is as high as 11.0°C at 3007 m a.s.l. (Han et al., 2008). The following vegetation 
types can be found around the glacier: alpine meadow (2900-3600 m a.s.l.) and alpine cushion 
plants (3600-4250 m a.s.l.) (Sabit et al., 2016). The land near the glacier terminus is barren and 
covered by massive glacial gravel, but an alpine meadow occurs at 2979 m a.s.l. (Fig. 2g—i). 


2.2 Unmanned aerial vehicle (UAV) survey and data collection 


The Mavic Pro (Fig. 3), which was employed in this study, is an electric rotary quadcopter 
manufactured by DJI Innovation Company, Shenzhen City, China. It is equipped with a 
high-precision vision positioning system that enables it to precisely capture images with geodetic 
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Fig. 2 Photographs of surface features on the Kekesayi glacier foreland (a, b, and c), Jiangmanjiaer glacier 
foreland (d, e, and f), and Koxkar Baxi glacier foreland (g, h, and i) along elevational gradients 


accuracy even in conditions having a poor global navigation satellite system signal. A three-band 
FC220 digital camera, integrated in the front of Mavic Pro, features a 12-megapixel resolution 
(4000x3000 pixels) and maintains a horizontal orientation when stationary. 

During the period from 26 September to 2 October, 2018, we selected 23 sampling sites along 
elevational gradients on the basis of the vegetation characteristics (FVC, plant height, and 
vegetation type) on the three glacier forelands. These sites had obvious differences in the plant 
community characteristics and could reflect the real-life feasibility of vegetation distribution. At 
each site, the UAV was taken from 2.0 m in height and covered the quadrats with 40 mx40 m in 
size (Fig. 3). Prior to each flight, waypoints were adjusted to follow the topography of the glacier 
forelands. The UAV took photographs vertically downward while flying autonomously between 
waypoints generated in the flight planning software. During the flight, the camera was triggered 
every 3 s as the UAV traveled at 3.0 m/s, and all images were stored in the JPEG format to aid 
later processing. The white balance of the camera was set in accordance with the actual weather 
conditions. At least 16 images were obtained on each flight under normal circumstances. The high 
spatial resolution (0.07 cm) of all images was beneficial for extracting the necessary underlying 
surface information. 

No meteorological station existed on the three glacier forelands, with the exception of one 
automatic weather station on the Koxkar Baxi glacier foreland. Hence, meteorological data 
(temperature and precipitation) were calculated through empirical formulas. On the Kekesayi and 
Jiangmanjiaer glacier forelands, the following equations were applied: 


A; — 701 
AP». sasore M, 


. 1 
amir 100 (1) 
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Fig. 3 Field experimental design of the sampling site selection and unmanned aerial vehicle (UAV) survey, as 
well as the later processing of aerial photographs. (a), distribution of sampling sites along elevational gradients; 
(b), photograph showing the UAV; (c), design of UAV survey; (d), aerial photograph; (e), interpretation of aerial 
photograph. 


AMan T E A 2) 
where APpamir and AMTpamir denote the annual precipitation (mm) and annual mean temperature 
(°C) of the eastern Pamir Plateau, respectively; values of 13.7 and 605.0 are the precipitation 
lapse rate (mm/100 m) and annual precipitation (mm) at 7010 m a.s.l., respectively (Duan et al., 
2007); values of 0.71 and 0.70 are the temperature lapse rate (°C/100 m) and annual mean 
temperature (°C) at 3250 m a.s.l., respectively (Luo, 1994; He et al., 2005); and A; represents the 
elevation variable (m a.s.1.). 


On the Koxkar Baxi glacier foreland, the equations used are as follows: 


AP rianshan = 069.4 + 49.8 x rey E 
100 
AMT rianshan = 0.70 — 0.63 x Aaa M 


where APrianshan and AMTrianshan reflect the annual precipitation (mm) and annual mean 
temperature (°C) of the western Tianshan Mountains, respectively; values of 49.8 and 669.4 are 
the precipitation lapse rate (mm/100 m) and annual precipitation (mm) at 3000 m a.s.l., 
respectively (Li et al., 2012); and values of 0.63 and 0.70 are the temperature lapse rate (°C/100 
m) and annual mean temperature (°C) at 2950 m a.s.l., respectively (measured data) (Han et al., 
2008). 

The calculated results are presented in Table 1. 
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Table 1 Annual mean temperature (AMT) and annual precipitation (AP) on the three glacier forelands 


Glacier foreland Elevation (m a.s.l.) AMT (°C) AP (mm) Aspect 
Kekesayi glacier foreland 3843—4004 —4.7 to -3.5 171.1-193.2 Northeast 
Jiangmanjiaer glacier foreland 3631-4241 —6.3 to —2.0 142.1-225.6 Southwest 

Koxkar Baxi glacier foreland 2979-3099 —0.2 to 0.5 658.9-718.7 South 


2.3 Evaluation of fractional vegetation cover (FVC) and gravel cover (GC) 


For the purpose of accurately obtaining ground-cover data such as FVC and GC, we removed 
low-quality images and adopted the traditional maximum likelihood classifier to process 
high-quality ones by utilizing the ENVI software (Harris Geospatial Solutions, Boulder, CO, 
USA). The aerial photographs were positioned by Pix4D software (Pix4D Company, Lausanne, 
Switzerland) prior to classification. The deformation of the images was negligible due to low 
flight height (2.0 m) and relatively flat aerial sites. Individual photographs were classified 
independently using a supervised classification-based maximum likelihood classification method. 
According to the purpose of this study, we divided photographs into three land cover types, which 
are vegetation, glacial gravel, and others (such as water body, barren land, and so on) (Fig. 3). We 
conducted the classification accuracy assessments of the resulting layers of aerial photographs by 
contrasting the sample class of the classified layer with the reference layer. The overall accuracy 
and Kappa coefficient were computed so as to evaluate the degree of the classification accuracy 
of the error matrix (Muzein, 2006). The overall accuracy is the sum of correctly classified values 
(diagonals) divided by the total number of randomly generated reference values of the error 
matrix (Lillesand and Kiefer, 2000). The Kappa coefficient, which determines the difference 
between the actual agreement of classified map and the chance agreement of random classifier 
compared with reference data, was calculated as follows: 
k k 
Ks NY ita = za) (5) 


N? > es x X,) 


‘ hd k f 
where Kna is the Kappa coefficient; N presents the total number of values; Xab ÍS the 


k : 
observed accuracy; and pae x x) is the chance accuracy. 


The results were tested with a Kappa coefficient of over 0.58, which indicated that data quality 
was acceptable. Finally, we calculated FVC and GC based on the classified results, and the 
equations are as follows: 


Pixcel number, 


FVC (%)= x 100%, (6) 


Pixcel number, 


Pixcel number 
——_——+£ x 100%, (7) 
Pixcel number, 


GC (%)= 
where FVC is the fractional vegetation cover (%); GC is the gravel cover (%); Pixel number, is 
the pixel number of vegetation; Pixel number, is the pixel number of glacial gravel; and Pixel 
number, is the total pixel number within each aerial image. 

In this study, we regarded the mean FVC and mean GC calculated by 16 aerial photographs as 


reference data at each site. The flowchart for the calculation of FVC and GC is shown in Figure 4. 
2.4 Calculations of plant family composition, dominant species, and species richness 


The plant families were identified by aerial photographs (see Wei et al. (2021) for details), and the 
identified results were tested by the ground-based sampling. A total of 12 vascular plant families 
were identified in this study (Table 2). Besides, the dominant species within each sampling site 
were confirmed by the occurrence frequency as well as the cover of each species in 16 aerial 
photographs. This indicated that the dominant species had the highest occurrence frequency and 
coverage. Species richness is a common vegetation index that reflects species diversity in a 
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Fig. 4 Flowchart for the calculation of fractional vegetation cover (FVC) and gravel cover (GC) 


Kappa coefficient 


community or biotope (Karen et al., 2004). In this study, based on the characteristic difference of 
a variety of species, we visually identified and counted the species number by the 16 images 
within each site. After that, the species richness was obtained using the equation below: 


D = Species number appeared in unit area, (8) 


where D refers to the total number of species at each sampling site, namely, the species number 
appearing in 16 aerial photographs. 


Table 2 Plant family composition at sampling sites of the three glacier forelands 


Glacier foreland Elevation (m a.s.1.) Plant family composition 
4004 Cyperaceae and Compositae 
3963 Cyperaceae, Compositae, and Gramineae 
Kekesayi glacier 3913 Cyperaceae, Compositae, Gramineae, and Rosaceae 
foreland 3890 Cyperaceae, Compositae, Gramineae, Rosaceae, Ephedraceae, and Fabaceae 
3873 Cyperaceae, Compositae, Gramineae, Rosaceae, Ephedraceae, and Fabaceae 
3843 Cyperaceae, Gramineae, Tamaricaceae, and Polygonaceae 
4241 Gramineae and Compositae 
4206 Gramineae, Compositae, Fabaceae, and Caryophyllaceae 
4185 Gramineae, Compositae, Fabaceae, and Caryophyllaceae 
4149 Gramineae, Compositae, Fabaceae, Rosaceae, and Crassulaceae 
Jiangmanjiaer 4063 Gramineae and Compositae 
glacier foreland 4010 Gramineae and Compositae 
3793 Gramineae, Compositae, and Chenopodiaceae 
3730 Gramineae and Cyperaceae 
3661 Gramineae and Cyperaceae 
3631 Gramineae and Cyperaceae 
3099 No 
3049 No 
3011 Rosaceae 
oe d 2992 Gramineae, Cyperaceae, Polygonaceae, Crassulaceae, and Asteraceae 
2988 Gramineae, Cyperaceae, Polygonaceae, Crassulaceae, and Asteraceae 
2986 Gramineae, Cyperaceae, Polygonaceae, Crassulaceae, and Asteraceae 


2979 Gramineae, Cyperaceae, Polygonaceae, Crassulaceae, and Asteraceae 
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2.5 Linear mixed-effects model 


The linear mixed-effects model was applied to assess the changes in FVC as a function of GC 
(Meng et al., 2007). Fixed and random effects, as the major components of the model, had to be 
determined prior to calculation. Specifically, we first analyzed the effect of GC on FVC; 
therefore, GC could be considered as the fixed effect. At the glacier basin scale, aspect affected 
soil hydrothermal conditions via altering solar irradiation and subsequently altered the vegetation 
pattern. At the regional scale, the regional climate was one of the key elements influencing 
vegetation development. This led to the identification of regional climate and aspect as random 
effects. The analyses were conducted using the Ime function in the nlme package in the statistical 
program R v. 3.6.1 (https://cran.r-project.org/). Finally, we used the r.squaredGLMM in the 
MuMIn package in the statistical program R v. 3.6.1 for the purpose of assessing the accuracy of 
the model. The function provided two measures: R’m and R’c. The first reports the R? of the 
model with just fixed effects, whereas the second represents the R? of the full model (including 
fixed and random effects). R? is the criteria indicating the fitness of the models. If R? is similar 
between models, it is most importantly that Rc is only slightly different compared to R?m, which 
implies that the inclusion of random effects does not improve the accuracy. 


3 Results 


3.1 FVC along elevational gradients 


On the Jiangmanjiaer and Kekesayi glacier forelands of the eastern Pamir Plateau, FVC varied in 
an inverted "N-shaped" pattern with increasing elevation, but it exhibited declining trends with 
increasing elevation in general (Fig. 5). FVC was highest on the bottom of valleys (old sites), 
whereas it was lowest near the glacier terminus (young sites). Particularly, on the Jiangmanjiaer 
glacier foreland, the maximum and minimum of FVC were 91.9% (at 3631 m a.s.l.) and 4.8% (at 
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Fig. 5 Variations in percentages of land cover types (glacial gravel, vegetation, and others) along elevational 
gradients and the relationships between GC and FVC on the Kekesayi glacier foreland (a, d, and g), Jiangmanjiaer 
glacier foreland (b, e, and h), and Koxkar Baxi glacier foreland (c, f, and i). VC, vegetation cover. 
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4241 m a.s.l.), respectively, whereas those on the Kekesayi glacier foreland were 65.7% (at 3843 m 
a.s.l.) and 35.2% (at 3913 m a.s.l.), respectively. Furthermore, the spatial variation of FVC on the 
Jiangmanjiaer glacier foreland was more dramatic in contrast to that on the Kekesayi glacier 
foreland, with a range of 87.1%, while the Kekesayi glacier foreland was densely vegetated (Fig. 5). 

On the Koxkar Baxi glacier foreland of the western Tianshan Mountains, FVC exhibited a 
linear decline trend with increasing elevation (Fig. 5c and f). Just like the spatial changes in FVC 
on the two glacier forelands of the eastern Pamir Plateau, its maximum (86.9%), as well as the 
minimum (0.0%), also occurred on the valley floor (old sites) and at the glacier terminus (young 
sites). FVC on the Koxkar Baxi glacier foreland of the western Tianshan Mountains showed a 
more pronounced growth trend with decreasing elevation compared with that on the two glacier 
forelands of the eastern Pamir Plateau (Fig. 5). 


3.2 Plant family composition along elevational gradients 


On the Kekesayi and Jiangmanjiaer glacier forelands, despite the fact that the surveyed glacier 
forelands were covered by cold-desert vegetation, the types and number of plant families 
displayed obvious differences along elevational gradients. On the Kekesayi glacier foreland, 
vegetation was primarily composed of individuals of the following families: Gramineae, 
Compositae, Cyperaceae, Ephedraceae, Fabaceae, Rosaceae, Tamaricaceae, and Polygonaceae. 
Cyperaceae was the dominant plant family and appeared at all sampling sites, whereas 
Tamaricaceae and Polygonaceae were the rarest plant families, only present near the bottom of 
valleys (old sites; Table 2). From the elevation of 4004 m a.s.l. to the glacier terminus (young 
sites), the land was covered by a lot of loose glacial gravel and sparse vegetation; thus, the 
simplest plant family composition appeared there. Moreover, the near valley bottom (3873 m 
a.s.l.) had the richest plant families, comprising 75.0% of all plant families. Plant family 
composition at the valley bottom (old sites) was simple and relatively stable; four plant families 
existed in this area, namely Gramineae, Cyperaceae, Tamaricaceae, and Polygonaceae. On the 
Jiangmanjiaer glacier foreland, the plant families consisted of Gramineae, Compositae, 
Cyperaceae, Fabaceae, Caryophyllaceae, Rosaceae, Crassulaceae, and Chenopodiaceae. 
Gramineae was the leading family and occurred at all sampling sites (Table 2). The variation in 
plant family composition along elevational gradients was stochastic. The richest plant family 
composition appeared at 4149 m a.s.l., whereas the simplest plant family composition occurred 
from the elevation of 4241 m a.s.l. to the glacier terminus (young sites). Even though each 
sampling site contained specific dominant species, Artemisia minor and Carex moorcroftii were 
the dominant species found at all sampling sites on the Kekesayi and Jiangmanjiaer glacier 
forelands, respectively (Fig. 6a and b). 

On the Koxkar Baxi glacier foreland, overall six plant families were recorded, including 
Rosaceae, Gramineae, Cyperaceae, Polygonaceae, Crassulaceae, and Asteraceae. From the 
elevation of 3049 m a.s.l. to the glacier terminus (young sites), there was no vegetation on the 
barren land; however, there was plenty of glacial gravel. Plants did not appear until the elevation 
of 3011 m a.s.l., in which only plants of Rosaceous family were distributed sporadically in the 
surroundings of glacial gravel or in crevices among them (Fig. 2h). With lowering elevation, plant 
family composition became increasingly abundant. Finally, a stable alpine meadow appeared on 
the valley floor (2979 m a.s.l.), where plant family composition was the richest and consisted of 
Gramineae, Cyperaceae, Polygonaceae, Crassulaceae, and Asteraceae (Fig. 2i; Table 2). Kobresia 
pygmaea and Carex melanantha were widely distributed below 2992 m a.s.l. (Fig. 6c). 


3.3 Species richness along elevational gradients 


On the Kekesayi and Koxkar Baxi glacier forelands, species richness first increased and 
subsequently decreased with increasing elevation (Fig. 7). Their maximum and minimum 
occurred near the bottom of valleys (old sites) and near the glacier terminus (young sites), 
respectively. On the Jiangmanjiaer glacier foreland, species richness displayed a fluctuant change 
of decreasing firstly, then increasing, and finally decreasing with rising elevation. The sites where 
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the extremum occurred were totally opposite to those on the other glacier forelands (Fig. 7). 

Overall, species richness exhibited decreasing trends on the Kekesayi and Koxkar Baxi glacier 
forelands with increasing elevation; however, it showed an increasing trend on the Jiangmanjiaer 
glacier foreland with rising elevation (Fig. 7). 


4 Discussion 


4.1 Community development of vascular plants on the three glacier forelands 


In the early stage of plant community development, vascular plants on the all glacier forelands 
shared a similar picture of low FVC and species diversity (Burga et al., 2010; Schumann et al., 
2016), whereas this landscape started to differ in the middle and later succession stages among the 
glacier forelands as a result of multiple factors, such as natural disturbances, terrain age, regional 
climate, topographic conditions, and so on (Tishkov, 1986; Jones and Roger, 2005; Raffl et al., 
2006; Dolezal et al., 2008; Burga et al., 2010; Wietrzyk et al., 2018). As a matter of fact, this 
natural scene also occurred on the three surveyed glacier forelands. In the newly exposed surfaces 
after the glacier retreated, certain microorganisms first colonized on the barren land, and 
eventually their dead organic matter facilitated the colonization of pioneer plants by supplying the 
initial soil as well as nutrient conditions. Accordingly, the early stage of plant community 
development had the lowest FVC and species abundance, which was in accordance with earlier 
studies on the Urumqi Glacier No. 1 foreland of northwestern China (Wei et al., 2021), 
Skaftafellsj6kull glacier foreland of southern Iceland (Glausen and Tanner, 2019), and glacier 
foreland in the European Alps (Schumann et al., 2016). These colonized pioneer plants further 
promoted the soil development and the accumulation of nutrients in soil (Wietrzyk et al., 2018), 
and thereafter vascular plants began to expand as succession proceeded. The spatial changes in 
vascular plant characteristics (FVC, plant family composition, and species diversity) had evident 
differences under different topographic (elevation and aspect) and climatic conditions. 
Specifically, despite the fact that FVC on the Kekesayi and Jiangmanjiaer glacier forelands 
located in the same glacier basin presented a similar variation trend along elevational gradients, 
the shady slope (Kekesayi) was more densely vegetated in contrast to the sunny slope 
(Jiangmanjiaer). Moreover, the Kekesayi glacier foreland had formed small vegetation patches, 
whereas vegetation on the Jiangmanjiaer glacier foreland tended to grow individually with the 
exception of the bottom of the valley where a grassland occurred. This was consistent with 
previous literature findings (Ostendorf and Reynolds, 1998; Pearson et al., 1999; Schumann et al., 
2016). Furthermore, plant species diversity on the Kekesayi glacier foreland was more abundant 
compared with that on the Jiangmanjiaer glacier foreland. This pointed out the fact that change in 
solar radiation caused by aspect might accelerate heat and moisture migration of topsoil, and 
subsequently could alter the spatial pattern of vegetation colonization and development. Our 
results were consistent with previous studies (Ostendorf and Reynolds, 1998; Schumann et al., 
2016), which indicated that aspect affects communities through altering solar insolation and hence 
temperature and effective moisture, thereby impacting biodiversity and composition. 

In addition to topographic factors, regional climate was another critical factor affecting 
vegetation dynamics. However, its effect on vegetation had usually been neglected except for the 
studies of Robbins and Matthews (2010) and Schumann et al. (2016). Meteorological data (Table 
1) revealed that climatic conditions on the survey glacier forelands of the eastern Pamir Plateau 
were colder and drier than those on the survey glacier foreland of the western Tianshan 
Mountains. Naturally, the difference in climatic conditions was also reflected in the ecosystem 
structure in the later stage of plant community development. In particular, on the relatively warm 
and humid glacier foreland (Koxkar Baxi), elevational zonality of vegetation was obvious, and an 
alpine meadow (at 2979 m a.s.l.) emerged just 576 m from the glacier terminus (Fig. 2i). On the 
other hand, the cold and dry glacier forelands of the eastern Pamir Plateau had sparse and 
scattered vascular plants and only progressed toward alpine grassland. This finding was consistent 
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with the study on the glacier forelands of the eastern and western Alps (Schumann et al., 2016). 
The study on the European Alps, nevertheless, contended that the differences in vegetation on the 
glacier forelands of the eastern and western Alps could not be explained by distinct climatic 
conditions, but rather by different species pools and treeline elevation (Schumann et al., 2016). 


4.2 Effects of glacial gravel on plant colonization and development on the three glacier 
forelands 


In the harsh environment of glacier forelands, vegetation colonization and development were 
primarily dominated by multiple biotic and abiotic factors, which made vegetation colonization in 
the early succession stage extremely difficult, and only a few pioneer plants had an opportunity to 
colonize in these barren lands (Dong et al., 2016). The three glacier forelands that were surveyed 
likewise experienced a similar situation. Interestingly, these pioneer species usually grew on the 
sides of glacial gravel or in the gaps, which was especially common on the newly exposed 
surfaces. Previous studies highlighted that this occurrence revealed how plants adapted to their 
harsh environment, but these studies only provided a qualitative account of how glacial gravel 
aided in the colonization and growth of vegetation (Thuiller et al., 2005; Mondoni et al., 2015). 

In this study, we quantitatively explored the function of glacial gravel in the colonization and 
growth of vascular plants on the glacier forelands. In the entire stages of plant community 
development, FVC was negatively related to GC (Fig. 5g—i; Table 3), but their relationships 
varied depending on the different topographic and climatic conditions. Specifically, FVC on the 
glacier forelands of the eastern Pamir Plateau and western Tianshan Mountains decreased in a 
power function and a linear function with increasing GC, respectively. In addition to the spatial 
relationships between vascular plants and glacial gravel (Figs. 2a and b, and 5g and h), we 
speculated that glacial gravel might facilitate early stage vegetation colonization and development 
by providing suitable hydrothermal conditions and preventing damage of herbivores, but it might 
also impose potential restrictions on vegetation expansion in the middle and late stages by 
reducing available space. Similar findings from earlier investigations have been reported 
(Stécklin and Bäumler, 1996; Jumpponen et al., 1999; Niederfriniger Schlag and Erschbamer, 
2000; Mong and Vetaas, 2006). According to Table 3, R’c was significantly different from R?m 
when we analyzed the changes in FVC as a function of GC on the glacier forelands of the eastern 
Pamir Plateau, suggesting that aspect could improve the model accuracy. This finding indicated 
that the effect of aspect on vegetation should not be neglected when analyzing the role of glacial 
gravel in vegetation development at the basin scale. In addition, it has been found that Rc was 
higher than R’m when evaluating the changes in VC as a function of GC in the eastern Pamir 
Plateau and western Tianshan Mountains (Table 3), implying that the effect of glacial gravel on 
vegetation distribution was mediated by regional climate at the regional scale. 


Table 3 Effects of gravel cover (GC) on fractional vegetation cover (FVC) at the glacier basin and regional 
scales using linear mixed-effects model and the assessment of the model accuracy 


Response Region AIC BIC logLik Intercept Slope Rk? Rm Rc 
Change in FVC as a EP 144.15 147.98 66.07 54.82 3.76* -0.53 0.21 0.59 
function of GC EP-WT 212.45 218.72 -100.23 63.30 6.15" -0.74 0.37 0.43 


Note: EP represents the eastern Pamir Plateau; EP-WT represents the eastern Pamir Plateau and western Tianshan Mountains; 
AIC, BIC, and logLik are the Akaike, Schwarz-Bayesian Information Criterion, and logarithmic likelihood used for model 
selection, respectively; Intercept and Slope are the two parameters of the linear mixed-effects model; R? is the criteria indicating 
the fitness of the models; R?m and R?c are the two measures of the assessment of the model accuracy. *, P<0.05 level. 


4.3 Advantages and limitations of UAV in the investigation of vegetation on the glacier 
forelands 


The UAV has been extensively employed in studying vegetation dynamics in recent years 
(Dunford et al., 2009; Chen et al., 2016; Mead and Arthur, 2020). However, its application for the 
investigation of vegetation on the glacier forelands is rarely reported with the exceptions of Eichel 
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et al. (2017) and Wei et al. (2021). In this study, we investigated the vegetation data (FVC, plant 
family composition, and species richness) on the three glacier forelands using an UAV, which 
demonstrated five benefits that are mentioned in Section 1. However, this advanced approach also 
had its limitations, including species identification in the laboratory requiring lots of time, low 
quality of aerial photographs affecting the outcomes of species identification, and specific plants 
(creeping and low-growing plants, and plants that matched ground colors) increasing difficulties 
in FVC extraction and species identification (Wei et al., 2021). Furthermore, the inherent 
advantages of UAV in monitoring vegetation were further constrained by adverse field time, 
which was also evident in this study. 


5 Conclusions 


This study aimed at assessing the spatial differences of vascular plant community development on 
the three glacier forelands (Kekesayi, Jiangmanjiaer, and Koxkar Baxi) under various topographic 
and climatic conditions, as well as exploring the effects of glacial gravel on plant community 
development by integrating the ground-based sampling and aerial photography. The findings 
pointed out that the community characteristics of vascular plants exhibited different spatial 
changes along elevational gradients at different spatial scales. Specifically, at the glacier basin 
scale, the shady slope (Kekesayi) was more vegetated in comparison to the sunny slope 
(Jiangmanjiaer); at the regional scale, the relatively warm and humid glacier foreland (Koxkar 
Baxi) had the highest FVC. Meanwhile, species diversity of vascular plants on the glacier 
forelands gradually decreased as elevational gradients increased, with the exception of that on the 
Jiangmanjiaer glacier foreland. Moreover, the correlations between GC and FVC showed that the 
role of glacial gravel in the colonization and development of vascular plants varied as community 
succession proceeded. Namely, glacial gravel tended to play a positive role in plant colonization 
of the early stage; however, it played a negative role in the plant community development of the 
middle and later stages. Numerous other factors, such as topography (elevation and aspect) and 
climate, also had an impact on this phenomenon. These findings provided a strong implication for 
ecological restoration in glacier areas with harsh environment. 
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